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I Abstract. In non-relativistic Brueckner calculations of nuclear matter, the self-consistent single 

CN ■ particle potential is strongly momentum dependent. To simplify the calculations, a parabolic 

^ , approximation is often used in the literature. The variation in the binding energy value intro- 

duced by the parabolic approximation is quantitatively analyzed in detail. It is found that the 
approximation can introduce an uncertainty of 1-2 MeV near the saturation density. 
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' 1. Introduction 

o ■ 

^ ' It is one of the fundamental issue in nuclear physics to evaluate the nuclear matter binding 

energy and saturation properties, starting from a realistic nucleon-nucleon (NN) interac- 
tion with no free parameter. This old project requires the solution of a complex many- 
body problem, and has received several contributions and improvements along the years, 
^ I beginning as far back as the middle of the last century. One of the main approaches to 

^ I this long standing problem is the so-called hole-expansion or Bethe-Brueckner-Goldstone 

(BEG) theory [|I|. The first real breakthrough in this scheme was the introduction of the 
^ . self-consistent single particle potential at the two hole-line level of approximation, which 

is then usually referred as the Brueckner-Hartree-Fock (BHF) approximation |I|. The 
introduction of the self-consistent potential drastically improves the results. In particular, 
the binding energy and saturation density, which otherwise would turn out unreasonable, 
move to values which can be considered an acceptable starting approximation. The re- 
maining discrepancy could be summarized in the celebrated "Coester band" [Q, along 
which the saturation points for different NN interactions were clustering and which miss- 
es the empirical region ( corresponding approximately to a binding energy per nucleon 
of -16 MeV and a nucleon density of 0.17 fm~^ ). Later, the Liege group stressed 
the relevance of the choice of the single particle potential. In particular they suggested 
the use of the "continuous choice", which indeed appears to move the saturation point 
towards the empirical one, but still missing it 0, ^. A period of major developments 
took place in the latest two decades. Starting from the works by B.D. Day 0], the hole- 
line expansion was analysed up to the three hole-line level of approximation. A strong 
indication of convergence of the expansion was obtained P, |^. Furthermore, BHF cal- 
culations with the continuous choice seem to get a substantially smaller corrections from 
three-body correlations 0. The results confirm that the empirical saturation point is 
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still missed, and therefore that three-body forces are needed in the nuclear hamiltonian 
[0]. In the meanwhile the relativistic Dirac-Brueckner (DB) method was developed [§], 
which already at the two hole-line level of approximation appears to be able to repro- 
duce the empirical saturation point. The main relativistic correction introduced by the 
DB method is due to the structure of the Dirac 4-spinors, which in the medium appear 
"rotated" with respect to the free ones. The non-relativistic three-body forces and this 
relativistic effect of the DB approach are probably two faces of the same dynamical effect 
[^]. The many-body theory has reached, therefore, such a precision that it is possible to 
test the nuclear hamiltonian. Because of that, the time seems to be appropriate to check 
the reliability of the approximations which are commonly employed in BBG calculations 
of nuclear matter. 

In this letter we consider the BHF in the continuous choice and we analyze quantita- 
tively the uncertainty of the results which comes out by approximating the single particle 
self-consistent potential with a parabolic form. This approximation is quite popular, since 
it allows to calculate the potential, at each iteration, only for few momenta, thus reducing 
drastically the computer time. The single particle potential, as obtained from fully self- 
consistent BHF calculations, is indeed strongly momentum dependent and not necessarily 
so simple parabola. 

2. Sketch of the formalism 

In the BHF approximation, the nuclear matter total energy E is obtained from the Brueck- 
ner G-matrix G{uj) according to the equation 

^= E ^ + ^ E {kih\Gie,,+e,,)\hh)A (1) 

ki<kp ki,k2<kp 

with \kik2)A = \k1k2) — \ k2k1) . Here fci? is the Fermi momentum, the summation over the 
momenta ki include spin and isospin variables. The single particle energies e^, appearing 
in the entry energy of the G-matrix, are given by 

e{k) = ^ + U{k) (2) 

where the single particle potential U{k) is determined by the self-consistent equation 

Uik)= {kk'\Giek,+ek,)\kk') (3) 

k'<kp 

The self-consistency is coupled with the integral equation for the G-matrix 

{kik2\G{Lu)\ksk4) = {kik2\v\ksk4,) + 

(4) 

+ Ek',k'Sk,k2\v\k',k',) ^'"''"iiy!:®"^'^^^ {k',k',\G{uj)\ksk,) 

where Qpik) = 1 for k < kp and is zero otherwise. The product Q{k,k') = (1 — 
QF{k)){l — Qpik')), appearing in the kernel of Eq. (^), enforces the scattered momenta 
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to lie outside the Fermi sphere, and it is commonly referred as the "Pauli operator" . The 
self-consistent set of equations are usually solved by an iteration procedure. The G-matrix 
can be expanded in partial waves, according to the classification of two-nucleon channels 
[0]. To avoid coupling between different two-body channels, the Pauli operator Q, as 
well as the two-body energies e^^ + e^^ in the denominator, are averaged over the angle 
between the relative momentum q = (fcg — k'^)/2 and the total momentum P = k'^ + k'^. 
Despite this approximation, which has been tested recently in ref. |T^, the numerical 
solution of the coupled equations is quite time consuming, since the single particle 

potential U (k) must be calculated in a wide range of momenta with a fine enough grid. If 
one assumes that the potential U{k), or equivalently the single particle energy e{k), has 
approximately a quadratic form 

n'^k^ 

e(A:)-eo + ^ (5) 

then one can calculate the potential, at each iteration step, in few points only and in- 
terpolate the obtained values with a parabola. The approximation of Eq. @ is usually 
called the effective mass approximation, since then the spectrum has the same shape as 
the free one but with an "effective mass" m*. 

In order to test this approximation, we have performed a set of BHF calculations fully 
self-consistently without any assumption about the potential shape, as well as by forcing 
the potential to a parabolic shape by means a fitting procedure, and then compared the 
results. 



3. Results and discussion 

The performed BHF calculations include all two-body channels up to total angular mo- 
mentum J = llh. In a set of calculations we adopted the Argonne vig potential as 
the NN interaction. This potential belong to a new generation of realistic NN potentials, 
with an improved fit of the scattering data, which give similar results and cluster closely 



together in the Coester band |]I2[. The self- consistent single particle potential U{k) was 
calculated up to the momentum cut-off kmax = which turns out to be large 

enough in the considered density range [|1^]. The potential, for the Fermi momentum 
kp = 1.4/m~^, is displayed in Fig. 1 (full circles). It is numerically calculated with a grid 
step of 0.1/m~^ from the G-matrix, Eq. (^, and inserted as the entry potential at each 
iteration step, until convergence is reached, i.e. the potential and the binding energy are 
stable under iteration with good accuracy. Stability within few KeV of the binding was 
systematically reached. The numerical method is described in ref. 

The quadratic approximation, at each step of the iteration procedure, is introduced by 
fitting the potential up to a certain maximum momentum kpjT- For definiteness, we have 
considered in detail two choices, namely kpij- = 2kp and kpix = kmax- At each iteration 
step, the potential U{k) coming directly from the G-matrix calculation is fitted with a 
parabola, which is then used as the entry potential for the next iteration. Convergence is 
reached when both potentials remain stable under this procedure. In this procedure one 
obtains, therefore, two potentials, one calculated from the G-matrix with the parabolic 
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input, and one from the parabolic fit to this potential. Of course, if the the potential 
coming directly from the G-matrix were indeed parabolic, the two potentials would closely 
agree. For the choice /cf/t = '^kp , in Fig. 1 the two potentials at convergence are 
displayed. In principle, one can calculate the nuclear matter binding energy from both 
potentials, but the result will be in general slightly different. As one can see, the fully 
self-consistent potential, obtained without any fitting procedure, as specified above, does 
not coincide with anyone of the two previous potentials, and these differences give a 
quantitative indication of the uncertainty introduced by the parabolic approximation. 
The corresponding saturation curves are reported in Fig. 2a. The parabolic potential 
produces a saturation curve in fair agreement with the one reported e.g. in ref. [|I0| 



Around saturation the parabolic approximation introduces a shift in the binding of 1-2 
MeV. The two choices for the potentials, discussed above, give different binding, since U {k) 
is not really parabolic, and the fitting procedure introduces necessarily an approximation. 
Another uncertainty is coming from the choice of kpiTi as can be seen in Fig. 2b, where 
the results for kpjT = k^ax are reported. In this case the discrepancy are larger for lower 
density, since then the potential U{k) becomes indeed flatter at momenta below kp. A 
more complete account of the dependence on the fitting range is reported in Fig. 3, where 
the binding at kp = 1.4/m~^ is reported as a function of kpip- 

In all cases the saturation curves appear distorted, and the saturation point shifted. 
Even if in some cases the saturation point seems to be "improved", this does not have 
any physical meaning, since, anyhow, it is mainly a spurious effect. 

Completely similar results are obtained with the "old" potential Argonne V14 . 



In conclusion, we have shown that the parabolic approximation for the single particle 
potential U{k) in the self-consistent Brueckner scheme introduces an uncertainty of 1-2 
MeV near the saturation density, and therefore it cannot be used in accurate calcula- 
tions. The full momentum dependence has to be retained, which prevents the use of a 
constant effective mass approximation. However, the uncertainty is not dramatic, and for 
approximate estimates of nuclear binding it can be useful. 
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Figure captions 



Fig. 1.- Single particle potential as function of momentum. The full circles indicate 
the results of the fully self- consistent calculation, where the potential is taken at each 
iteration step as calculated from the Brueckner G-matrix. The solid line is the result of 
the parabolic approximation. The parabolic potential, used as input for the G-matrix, 
produces the potential indicated by the squares. 

Fig. 2.- Saturation curves as obtained for the fully s elf- consistent calculaiion (full circles), 
from the parabolic fit (triangles) and from the single particle potential obtained from the 
parabolic input (squares). Panel a corresponds to kpir — 2A;f, panel b to kpir — kmax- 

Fig. 3.- Dependence of the binding energy from the value of kpiT-i using either the parabolic 
fit potential (triangles), or the parabolic input (squares) potential. The arrows indicate 
the value of the binding energy obtained in the fully self- consistent calculation, where no 
fitting procedure is introduced. 
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